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We propose a theoretical framework within which a robust mechanical definition of precursors to 
plastic instabilities, often termed ‘soft-spots’, naturally emerges. They are shown to be collective 
displacements (modes) z that correspond to local minima of a ‘barrier function’ b(z). The latter is 
derived from the cubic expansion of the variation 5Us(s) of the potential energy upon displacing 
particles a distance s along z. We show that modes corresponding to low-lying minima of b(z) lead 
to transitions over energy barriers in the glass, and are therefore associated with highly asymmetric 
variations 5Uz(s) with s. We further demonstrate how a heuristic search for local minima of b(z) can 
a-priori detect the locus and geometry of imminent plastic instabilities with remarkable accuracy, 
at strains as large as 7 C — 7 ~ 1CU 2 away from the instability strain y c . Our findings suggest that 
the a-priori detection of the soft-spots field in model glasses can be effectively carried out by the 
investigation of the landscape of b(z). 


I. INTRODUCTION 


Plastic flow of disordered solids subjected to external 
loading is known to occur via localized rearrangements of 
small sets of particles, coined shear-transformations [T]. 
Such rearrangements have been identified in experiments 
on bubble rafts [2], foams [3], emulsions HUE and col¬ 
loidal glasses mi, as well as in atomistic computer sim¬ 
ulations of model glasses 7, 8, An example of such a 
shear-transformation, observed in a model glass in two 
dimensions deformed under athermal, quasi-static (AQS) 
shear jjjjj, is displayed in Fig. |TJo. Shear-transformations 
are known to self-organize in spatially correlated pat¬ 
terns [IMS] in solids subjected to large stresses and low 
deformation rates. Their densities and other statistical 
properties, and mechanical consequences, are a subject 
of much recent debate [ITIIM] . Two questions, central 
to theoretical descriptions of elasto-plasticity, that we 
address in this work, are: can shear-transformations be 
predicted a-priori, and, if so, how? 

The micro-mechanical process in which an athermal 
disordered solid destabilizes under quasi-static deforma¬ 
tion is understood, asymptotically close to an instabil¬ 
ity strain 7 C , as a saddle-node bifurcation of the po¬ 
tential energy U P3 I2Z3- The immediate precursors 
to shear-transformations at strains 7 —t 7 C are identi¬ 
fied as destabilized eigenfunctions T c (i.e. their associ¬ 
ated eigenvalues vanish at 7 C ) of the dynamical matrix 
M.ij = , where 27 denotes the coordinate vector 

of the i th particle. Such an eigenfunction is presented in 
Fig. 0. In the following, we refer to such eigenfunc¬ 
tions as destabilized modes , to distinguish them from the 
post-instability displacements of particles - agglomera¬ 
tions of shear-transformations - which can be spatially 
extended. In Fig. [I]) we demonstrate that, when the 
post-instability displacements are not spatially extended, 
but rather form an isolated elementary shear transforma¬ 
tion, their spatial structure is very similar to that of the 
destabilized mode. In contrast with the post-instability 
displacements that depend in general on a specific choice 



FIG. 1. (color online). A plastic instability in a sheared two- 
dimensional model glass, a) The destabilized eigenfunction 
b) An elementary shear-transformation: the post-instability dis¬ 
placements that followed the instability of panel a), c) Nonaffine 
displacements v calculated at £7 = 7 C — 7 = 0.007. This delocal¬ 
ized field is used as the initial conditions z- in i for the minimization 
of b(z) (see main text), the result of which is the plastic mode n 
displayed in panel f). Panels d) and e) are intermediate states 
along the minimization of b(z). 
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of dynamics El. and on external control parameters such 
as temperature [15] . strain tm, and strain rate |16j . the 
spatial stricture of destabilized modes is an intrinsic 
characteristic of the multi-dimensional potential energy 
function, and is therefore the focus of the present study. 

A robust mechanical definition of the precursors of 
plastic instabilities away from instability strains has not 
yet been put forward. Much effort has been dedicated re¬ 
cently to studying the role played by low-frequency nor¬ 
mal modes in determining these precursors P5b32j . One 
key difficulty encountered in such studies is that low- 
frequency plane-waves, which have no appreciable effect 
on plasticity 113. dominate the lower parts of the spectra 
of conventional model glasses, thus hindering attempts to 
use low-frequency modes to define flow-defect densities 
and correlate them with rates of plastic flow. 

Another difficulty, which has been largely overlooked 
in the context of elasto-plasticity, is that mere frequen¬ 
cies of normal modes are not indicative of their relevance 
to plastic processes. In fact, modes which lead to me¬ 
chanical instabilities (i.e. take the system over energy 
barriers and into neighboring inherent states) appear as 
eigenfunctions of the dynamical matrix only very close to 
plastic instabilities, giving rise to difficulties in their de¬ 
tection and statistical quantification. Here we show that 
the effective detection of such modes away from plastic 
instabilities necessitates the consideration of the degree 
of asymmetry associated to the variation SU of the po¬ 
tential energy upon displacing the particles along those 
modes. We provide a theoretical framework that nat¬ 
urally embeds a mechanical definition of the precursors 
to plastic instabilities, and which effectively accounts for 
the said asymmetry. 


II. THEORETICAL FRAMEWORK 


We begin the discussion by considering an athermal 
elastic solid, of N particles in d dimensions, and let z de¬ 
note a Nd dimensional unit vector, i.e. Zj ■ Zj = 1. Here 
and in what follows, repeated indices, labeling particles, 
are understood to be summed over, unless indicated oth¬ 
erwise. The coordinates x are displaced in the direction 
defined by z according to 5x = sz, and we expand the 
potential energy U as: 


SU g (s) = U~(s) -U 0 ~ \k;s 2 + ±r,s 3 , (1) 

where Uq is the energy of the minimum in which the sys¬ 
tem resides, n~ = Mij : Zjij is the stiffness associated 

A r\ 3 T 7 • /\ y\ • • 

to z, and Tg = ds .gg dSk '■ ZiZjZk is referred to in the fol¬ 
lowing as the asymmetry associated to z. Within this 
cubic expansion, stationary points occur at s = 0 and 
s*(£) = — 2 ^; s = 0 corresponds to the minimum in 
which the system resides, while s* represents the saddle 
point (energy barrier) that separates this minimum and 
a neighboring inherent state. We thus define the energy 


difference between these stationary points, within the cu¬ 
bic expansion, as our barrier function: 


b(z) = \k%s1 + gTiS 3 



( 2 ) 


We emphasize that b(z) is defined for a particular con¬ 
figuration of an elastic solid in mechanical equilibrium, 
and is a function of the direction z. It has a rough land¬ 
scape |33j : in this work we focus on directions if that 
correspond to local minima of b{z). We refer to these 
directions in what follows as plastic modes. From the 
definition of b(z) it is clear that plastic modes ff are as¬ 
sociated with small stiffnesses k* and large asymmetries 
t*; they can be found numerically by minimizing b(z) 
over directions z, starting from some initial direction z ini , 
as demonstrated in panels c)-f) of Fig. [I] Small b(z)’s 
should appropriately describe low saddle points (barri¬ 
ers) that separate the system from neighboring inherent 
states. We therefore expect modes ff that correspond to 
low-lying minima of b(z) (which are found by chosing an 
appropriate z iTli for the minimization), to encode infor¬ 
mation about imminent plastic instabilities. 




FIG. 2. (color online), a) Variations <5t/*(s) of the potential en¬ 
ergy upon displacing the particles a distance s along plastic modes 
it, obtained as described in the text. The curves correspond to 
5 7 = 10 -5 ,10 -4 , 3xl0 -4 ,8xl0~ 4 , 2xl0 -3 , 7xl0 -3 . b) Stiffnesses 
Kz = Mij '■ ZiZj associated to (circles) the plastic modes n used to 
calculate the variations plotted in panel a), and to (squares) the 
destabilized eigenfunction 'I' c , vs. S'y. 


A. Numerical demonstration 


The approach described above is demonstrated in 
Fig. [l] in panel a) we display a destabilized mode T c 
calculated at the firstly-encountered plastic instability in 
an athermally sheared model glass, here at a strain y c = 
0.011521. Prior to this instability, at strains 7 = 7 c — <5 7 , 
the nonaffine displacement field Vj = —Ad )) 1 • is cal¬ 
culated [27] . An example of v, calculated at S'y = 0.007, 
is shown in panel c). At this distance (in strain) from 
the instability, the nonaffine displacements v are largely 
delocalized. We use the normalized v = 'c/||f?|| as the 
initial conditions z ini for the minimization of b(z)\ snap¬ 
shots along the minimization are displayed in panels d) 
and e). Upon convergence, we find a local minimum 
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in the direction if, which is displayed in panel f). The 
resemblance between 7 r and the destabilized mode T c is 
striking: both the geometry and the core location appear 
to agree perfectly. 

This protocol is carried out over a broad range of in¬ 
tervals Sj, as specified in the caption of Fig. [5] For 
each £ 7 , after finding ft as described above, we calcu¬ 
lated its associated energy variation SUjt(s) and stiffness 
k*, which are displayed in panels a) and b) of Fig. [ 2 ] re¬ 
spectively. In this example, already at a distance of the 
order Sj ~ 10“ 3 to the instability strain, following the 
plastic mode it would carry the system above an energy 
barrier and into a neighboring minimum. 


III. DESTABILIZATION OF PLASTIC MODES 


We also plot in Fig. j2jo the stiffness k^, associated to 

the destabilizing mode Tc. We find that only very close 
to the instability (<5y < 10“ 5 ), the scaling ~ y/Sj 
holds [34], whereas the stiffness associated to if follows 
Kn~V&y up to strains of order 1 % away from the insta¬ 
bility. This finding supports the robustness of our defi¬ 
nition of plastic modes, and the usefulness of our frame¬ 
work. It also supports the picture proposed by a number 
of recent studies [HHM], that assumes the (reversible) 
destabilization process of a ‘soft spot’ in a deformed glass 
is predominantly coupled to the external load, and not 
to other coexisting (reversible) destabilization processes. 

The scaling can be derived as follows; modes 

it pertain to local minima of b(z) and therefore satisfy 
§1^=0, which implies that (see Appendix): 


d 3 U 


dxidx jdxk 


■ . „ T* . . 

: 7T,'7T fe = — Mi- 


"J ■ 


(3) 


Using this relation, we calculate the leading order varia¬ 
tion of the stiffness k * with strain as 


d,Kn 

dr) 


dMi i 

dj 


d 3 U 


7T i-Ma-M 


dxidxidxk 
d 2 U 


-1 


• 7T i^j^k 


dx k dj dSid 7 


d 2 U 


■(4) 


a2 t t 

As 7 —> 7c, k* -A 0, but 7777 • gg q goes to a constant, 
yielding the differential scaling relation ^ 1 ~ — -4-, and 
thus the observed scaling ~ y/5 7. 

The resolution of the plastic mode as seen in Fig. [l] 
uses the nonaffine displacements v as the heuristic guess 
for £ ini ; this choice is made to demonstrate the usefulness 
of the framework - despite the extended character of v, 
it has a large overlap with the plastic mode it, and thus 
resides in the basin of it on the landscape of b(z). Ob¬ 
taining the full field of plastic modes, however, requires 
using other heuristic z ini ’s, that reside in basins that be¬ 
long to other plastic modes. We leave the investigation 


of the optimal heuristics for the detection of the full field 
of plastic modes for future work. 



FIG. 3. (color online), a) Scatter plot of the barrier function 
Eq. § evaluated for eigenfunctions of A4, vs. their eigenvalues 
cj 2 . The eigenfunction ^ m i n represented by the circled data point 
is plotted in panel d), and used as the initial conditions £i n i for 
the minimization of the barrier function b(z); the resulting plastic 
mode 7 f is displayed in panel e). b) Variations 5Uz(s ), calculated by 
displacing the particles according to Sx = s^ m i n (continuous curve) 
and by Sx = sn (dashed curve), c) The products iV|T^|, averaged 
over bins of the participation ratio, see text for definitions. 


IV. COMPARISON TO NORMAL MODES 


How indicative are normal modes of imminent plastic 
instabilities, compared to plastic modes? In Fig. [3jx we 
present a scatter-plot of the barrier function evaluated 
at normal modes vs. the square of their associated 
frequencies w 2 , calculated for a few tens of undeformed 
(isotropic) solid realizations. A clear trend appears: 
smaller values of 6(4' [i ,) are found for lower-frequency 
modes. The circled data point represents the mode 'k min 
associated to the lowest value of 6 (’F LJ ) amongst all modes 
calculated; it is displayed in panel d). Remarkably, this 
normal mode displays the same spatial features as ob¬ 
served for destabilized modes, reinforcing that b(z) is 
indeed sensitive to ‘plastic-like’ modes. The variation 
SUy . (s) is plotted in panel b) (continuous line). De¬ 
spite possessing the smallest b amongst our entire ensem¬ 
ble of modes, 5U^, . (s) displays only a slight asymmetry 
between positive and negative displacements s, and the 
energy monotonically increases with |s|. Using 4 , min as 
the initial condition £ ini for the minimization of b(z), we 
find the plastic mode it displayed in panel e). On the face 
of it, T min and 7 r appear to be very similar in their spatial 
structure and geometry. However, examining the corre¬ 
sponding variation 6U#(s), represented by the dashed line 
in panel b), reveals a dramatic difference between them: 
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following 7r takes the system over a energy barrier, to a 
neighboring minimum. 

We further utilize our ensemble of normal modes to 
study the relation between the degree of localization of 
modes and their associated asymmetries Tq,. A simi¬ 
lar analysis was carried out in [32] in the context of 
the unjamming point [35)457] . We quantify the degree 
of localization of a mode T via its participation ratio 
e = [A r JA('iq • ’lb) 2 ] -1 ; localized modes have e ~ iV -1 , 
whereas maximally delocalized modes have e ~ 1. In 
Fig. [3J: we plot the means [35], averaged over modes 
T with similar participation ratios, for systems of N = 
1024 and N = 4096. We find that for p artic ipation ra¬ 
tios e < 10 -1 , the asymmetries follow |t\j,| ~ (eN) -1 . 
This can be explained with the following simple model: 
if there are effectively N a non-zero components in a 
normal mode (0 < a < 1), normalization then requires 
that a characteristic non-zero component is of magni¬ 
tude H^iH ~ N~^. The participation ratio is then ex¬ 
pected to follow e ~ 4V a_1 (due to summing over posi¬ 
tive terms). Since the pairwise potential is short ranged, 
and the tensor elements Q-F a - are of either sign, 
then Tq, consist of a sum over N a terms, each of order 
1111 3 ~ N~ ^, of random signs, and we therefore expect 
rq, ~ N~ a ~ (e-/V) _1 , in consistency with our measure¬ 
ment. For participation ratios e > 10 _ , this relation 
breaks down, and asymmetries are much smaller than 
what is predicted by this simple model, which assumes 
that normal modes are random objects. Nevertheless, the 
same trend remains unchanged: delocalized modes are 
associated, on average, with more symmetric variations 
of the energy. These observations explain the localized 
nature of plastic instabilities found in deformed glasses, 
as can be seen, e.g., in Fig. Ek 



FIG. 4. (color online), a) Spatial decay of plastic modes, see text 
for definition. All modes analyzed decay as r 1— , as indicated by 
the continuous lines, where r is the distance to the core center. b) 
A plastic mode found by chosing a random z i n i. c) A plastic mode 
calculated in a disordered network of relaxed Hookean springs. d) 
A plastic mode found in a Lennard-Jones glass under isotropic ten¬ 
sion. 


random z ini . We find that this mode shares the same geo¬ 
metric features as the destabilized inodes 4 , c upon shear- 
induced plastic instabilities - a disordered core, and a 
long-ranged affine quadrupolar shear-like displacement 
field away from the core H2HZ3EZI- We thus conclude 
that plastic modes fr associated to different local minima 
of b(z) share similar structrural features, that do not de¬ 
pend on the particular minima to which they correspond. 


A. Effects of loading conditions 


V. STRUCTURE OF PLASTIC MODES 

To characterize the spatial structure of plastic modes, 
we define z 2 (r) as the median of the squared magnitude of 
the components z-i■ z, t (no summation implied), taken over 
a shell of thickness on the order of the nearest-neighbor 
distance, and of radius r away from the core of the plastic 
mode (detecting the locus of the core is explained in the 
Appendix). In Fig. |4^, we compare the spatial decay of 
two plastic modes, one obtained by setting 2 ini to be the 
direction of the nonaffine displacements (see definition 
above and Fig. [l]), and the other by setting z ini to be a 
random direction. These decay profiles are also compared 
to that of a destabilized mode '& c . We also show the 
decay profile of a plastic mode calculate in a 3D solid. 
We find that at distances r away from the core, plastic 
modes decay as r 1 ~ d . Remarkably, this is the same decay 
law found for the linear responses of displacements to 
dipolar point forces [35, 3])). 

In Fig. [4 ]d we present a plastic mode obtained with a 


We finally examine how the geometry of plastic modes 
depends on the loading conditions imposed on the solid. 
In panels c) and d) of Fig. [4j two additional examples 
of plastic modes fr obtained from a random z iai are dis¬ 
played; ft of panel c) was calculated in a disordered net¬ 
work of relaxed Hookian springs (all springs are neither 
stretched nor compressed) with an average of 4.1 springs 
connected to each node. It displays a similar spatial 
structure as that of plastic modes found in model glasses 
that are prestressed, i.e. in which finite forces are ex¬ 
erted between the constituent particles [40]. Our find¬ 
ings indicate that proximity to prestress-induced micro- 
mechanical buckling instabilities m is not the origin of 
the generic structure of plastic modes. 

The plastic mode if of panel d) of Fig. [4] was calcu¬ 
lated in a Lennard-Jones glass (with a pairwise poten¬ 
tial that includes an attractive term, see Appendix for 
details) under isotropic tension, just before macroscopic 
failure (here —p/B « 10~ 2 is at least 80% of the yield 
strain, where p is the pressure and B is the bulk mod- 
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ulus). We find in this case that in addition to the clear 
shear-like displacements that are typically seen in plas¬ 
tic modes found in glasses under compressive stresses, 
the dilatant part of the displacements due to the ten¬ 
sile loading conditions is apparent. We conclude that the 
loading conditions imposed on a solid can be reflected in 
the geometric features of its plastic modes; we leave the 
systematic study of this dependence for future work. 


VI. SUMMARY AND DISCUSSION 


In this work we demonstrated that modes 7r corre¬ 
sponding to local minima of the barrier function b(z), 
coined plastic modes, are indicative of directions in con¬ 
figuration space that lead to plastic instabilities, and 
more so compared to the most localized low-frequency 
normal modes. As such, our approach can serve as a 
solid basis for instability-detection algorithms. Such al¬ 
gorithms are highly desirable, as they can put to test 
theoretical frameworks of elasto-plasticity that involve 
the dynamics of a population of ‘soft-spots’. These al¬ 
gorithms need not be restricted to the investigation of 
plastic flow in disordered solids; the generality of our 
framework would render them suitable for studying a 
diverse set of systems, including dislocated crystalline 
solids, deeply supercooled liquids and proteins. 

Furthermore, our theoretical framework explains the 
origin of the localized nature of plastic instabilities. 
Building on our framework, we predict that the stiff¬ 
ness associated to plastic modes follows k ~ i/7 c — 7, 
and show numerically that this scaling holds over a large 
range of strains away from an instability strain 7 C . This 
adds relevance to recently proposed models that assume 
reversible destabilization processes of soft spots are de¬ 
coupled from each other. Finally, we have investigated 
the spatial features of plastic modes, and provide evi¬ 
dence that the detailed geometry of plastic modes is sen¬ 
sitive to the loading conditions imposed on the solid. 
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irnize the barrier function: 


^>4 


0 Z U ■ $.$. 

dxidx, ' 


d^u ' 

dxidxjdxk ' i' i > 


(Al) 


over directions z in the Ad-dimensional configuration 
space. Here U(x) is the potential energy, 27 denotes the 
Cartesian coordinate vector of the i th particle, and re¬ 
peated indices are understood to be summed over. Min¬ 
imizations are carried out using a standard nonlinear 
conjugate gradient minimization algorithm, which was 
modified such that the norm ||2|| remains of order unity, 
to maintain numerical stability. Notice the invariance 
b(z) = b(sz) for any finite length s, which allows us to 
control the norm without effecting the minimization pro¬ 
cess. 

Panels c)-f) of Fig. 1 in the main text show an exam¬ 
ple of a minimization of the barrier function. Panel c) 
is the initial guess z iai = v where v denotes the nonaffine 
displacement field, as defined in the main text. Panels d) 
and e) displays states encountered along the minimiza¬ 
tion algorithm, and panel f) is the plastic mode obtained 
upon convergence of the minimization algorithm to a lo¬ 
cal minimum of b(z). 


2. Mode spatial structure analysis 


The analyses of the spatial decay of modes £ were car¬ 
ried out as follows. First, we find the center of the core 
center , which we determine following 


x, 




E ae nl4ll 2 ^ 

EienPill 2 


(A2) 


where Q denotes the set of the w particles with the high¬ 
est ||zjj| 2 . The choice w = 4 was found to be the optimal 
choice for most of the plastic modes we analyzed, and we 
never used w > 10. We then divided the area of the sys¬ 
tem into rings (or shells in the 3D case) of radius r away 
from the center. The thickness of the rings was set such 
that there were at least 10 particles in each ring. For 
each ring, we found z 2 (r) defined as the median value 
over all ||ij|| 2 in that ring. 


Appendix A: Models and numerical methods 

1. Calculation of plastic modes 

We begin with describing our approach to calculating 
plastic modes; starting from an initial guess z ini , we min- 


3. Model definitions 

We continue with describing the models used in this 
work, starting with the model glass for which most data 
are presented: a 50:50 binary mixture of ‘large’ and 
‘small’ particles of equal mass to, interacting via radially- 
symmetric purely repulsive inverse power-law pairwise 
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potentials, that follow 


V^ipl (r ij ) 


£ 




(A3) 


where r t j is the distance between the i th and j th parti¬ 
cles, e is an energy scale, and x c is the dimensionless 
distance for which <^ IPL vanishes continuously up to q 
derivatives. Distances are measured in terms of the in¬ 
teraction lengthscale A between two ‘small’ particles, and 
the rest are chosen to be Ay = 1.18A for one ‘small’ and 
one ‘large’ particle, and Ay = 1.4A for two ‘large’ parti¬ 
cles. The coefficients C 21 are given by 


C2l 


(- 1 Y +1 (n + 2q)U {n+n) 

(2q — 2t)\\(2t)\\ (n — 2 )!!(n + 2£) c 


(A4) 


We chose the parameters x c = 1.48, n = 10, and q — 3. 
The density was set to be N/V = 0.86A -2 for our 2D 
systems, and N/V = 0.82A~ 3 for our 3D systems. We 
note that with these densities two coordination shells 


fall within the interaction range. Both the 2D and 3D 
systems undergo a computer-glass-transition at temper¬ 
atures of about T g ss 0.5 e/ks- Solids were created by 
a fast quench from the melt to a target temperature 
T <C T g , followed by an energy minimization using a 
standard nonlinear conjugate gradient algorithm. We 
employed two-dimensional systems of N = 4096 for pro¬ 
ducing the data presented in figures 1, 2, 3, and 4b of 
the main text. We measured decay profiles of plastic 
modes (plotted in Fig. 4a of the main text) in 2D sys¬ 
tems of N = 409600, and in 3D systems of N = 4096000. 
Normal mode analyses were carried out using the linear 
algebra package LAPACK [32] . and the numerical anal¬ 
ysis software MATLAB [43]. Systems were deformed by 
imposing a simple shear, namely each particle was dis¬ 
placed according to 


Xi -A Xi + S'yyt , (A5) 

2 (A 6 ) 

where Xi , yi denote the x and y coordinates of the i th 
particle. The strain increments 6 7 used were not larger 
than 10 -5 . 


Plastic modes were also calculated in two more systems: 2D disordered networks of relaxed Hookean springs (see 
panel (a) of Fig. [5] in this SM), which were created as described in [39] . and in 2D, 50:50 binary mixtures of particles 
interacting via the Lennard-Jones potential (see Fig. 4 of main text): 






+c4 te) 

0 


C 2 



121 . < qn 

\ • • — 

-W 1 


jT> x c 


(A7) 


Here x c = 2.5 is the dimensionless distance for which ip LJ vanishes continuously up to two derivatives, which implies 
that C 4 = 0.0006201261686784, c 2 = —0.00970155098112, and c 0 = 0.040490237952. The interaction lengths Ay- 
are identical to those chosen for the inverse-power-law system described above, as is the choice of units of length A, 
and e is an energy scale. Systems were first created at the initial density of N/V = 0.615A” 2 , which resulted in a 
small but positive pressure of p ss 0.05e/A 2 . We then applied small expansive strain increments <$ 7 , by applying the 
transformation 27 —> (1 + <$ 7)27 to the coordinates 27 , and subsequently minimizing the energy. We find that the 
system fails via cavitation at tensile pressures of p fail « —0.22e/A 2 . The bulk modulus [33] in the pre-failure states 
is B ss 18e/A 2 , which translates to a yield strain of p tail /B « 10 -2 . An example of a pre-failure solid is displayed in 
panel (b) of Fig. [5] in this SM. 


Appendix B: Shear coupling of plastic modes’ 
stiffnesses 


In this final section we provide an outline of the for¬ 
malism leading to Eq. (4) of the main text, for A 

detailed derivation will be presented elsewhere [45] . We 
begin with noting that, by construction, the barrier func¬ 
tion is invariant to inflations b(z) = b(sz), for any finite 
length s. This allows to define b(z) to be a function of 
the independent variables which is the form assumed 
below. 

r\ 2 j T 

For the sake of brevity we denote Mtj = d ~.g S . and 


Tjm — a-'u 
^ijk dxidxjdxk ’ 


and write the gradient of b(z) as 


db 

dzi 


= 4- 


M t 


- *u!" 


ijk : z j z k 


(Bl) 


Modes 7 r that correspond to local minima of b satisfy 
= 0, which immediately yields Eq. (3) of the 
main text: 


U% k : fijitk = / M, 

r 


(B2) 


We next turn to calculating the total derivative with re¬ 
spect to strain 7 of the stiffness k* as 


dj 


dMij 

d"f 


. „ „ . . a dir,- 

: iTiiTj + 2Mij : 7r *“^- 


(B3) 
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An equation for ^ is obtained by requiring that 7r re¬ 
mains a local minimum of b under the deformation, 
namely 


d 
d 'y 


db 

.dzi 


= 0. 


(B4) 


Using that 



d 2 3 4 5 6 7 8 9 10 11 b 

d 2 b 

d 2 b 


dzidzj 

nj dSidzj 

*=* ^ ^dZi 


■ TTi = 0, 

( B5 ) 

one finds that (i) ^ • n = 0, and that (ii) ||^|| goes to 
a constant at the instability strain 7 C . However, ^4 is 
singular, and we are thus left with 



dK* 

g?7 


dM.i 

dj 


TTiTTj 


(B6) 


as seen in Eq. (4) of the main text. Finally, following jJJ 
and references within, the total derivative with respect 
to strain of the dynamical matrix is 


dMij _ dMij + dMij dxk 

dj c >7 dx dj ’ 


(B7) 


where ^ = — A4 ij - 1 • are the nonaffine displace¬ 

ments, denoted v in the main text. Since as 7 C is ap¬ 
proached, ||u| I ~ (y c — 7) _5 as, to leading order we 
find 


FIG. 5. (color online), (a) An example of a two-dimensional disor¬ 
dered network of relaxed Hookean springs with a mean connectivity 
(z) £3 4.1, as used in this work, (b) An example of a pre-failure 
two-dimensional Lennard-Jones solid. 


dMij 

dj 


ijk ‘ V k 


as seen in Eq. (4) of the main text. 


(B8) 


[1] A. Argon, Acta Metall. 27, 47 (1979). 

[2] A. S. Argon and H. Kuo, Mater. Sci. Eng. 39, 101 (1979). 

[3] M. Dennin, Phys. Rev. E 70, 041406 (2004). 

[4] P. Hebraud, F. Lequeux, J. P. Munch, and D. J. Pine, 
Phys. Rev. Lett. 78, 4657 (1997). 

[5] J. Clara-Rahola, T. A. Brzinski, D. Semwogerere, 
K. Feitosa, J. C. Crocker, J. Sato, V. Breedveld, and 
E. R. Weeks, Phys. Rev. E 91, 010301(R) (2015). 

[6] V. Chikkadi, G. Wegdam, D. Bonn, B. Nienhuis, and 
P. Schall, Phys. Rev. Lett. 107, 198303 (2011). 

[7] D. Deng, A. S. Argon, and S. Yip, Proc. R. Soc. London, 
Ser. A 329, 549 (1989). 

[8] M. L. Falk and J. S. Langer, Phys. Rev. E 57, 7192 
(1998). 

[9] See Appendix for details about the models, units, and 
numerical methods used in this work. 

[10] C. Maloney and A. Lemaitre, Phys. Rev. Lett. 93, 016001 
(2004). 

[11] A. Tanguy, F. Leonforte, and J.-L. Barrat, 


Eur. Phys. J. E 20, 355 (2006). 

[12] C. E. Maloney and A. Lemaitre, Phys. Rev. E 74, 016118 
(2006). 

[13] S. Karmakar, E. Lerner, and I. Procaccia, Phys. Rev. E 
82, 055103(R) (2010), 

[14] K. M. Salerno, C. E. Maloney, and M. O. Robbins, 
Phys. Rev. Lett. 109, 105703 (2012). 

[15] J. Chattoraj and A. Lemaitre, Phys. Rev. Lett. Ill, 
066001 (2013). 

[16] A. Lemaitre and C. Caroli, Phys. Rev. Lett. 103, 065501 
(2009). 

[17] H.G.E. Hentschel, S. Karmakar, E. Lerner, and I. Pro¬ 
caccia, Phys. Rev. E 83, 061101 (2011). 

]18] E. Bouchbinder and J. S. Langer, Phys. Rev. Lett. 106, 
148301 (2011). 

[19] J.-C. Baret, D. Vandembroucq, and S. Roux, 
Phys. Rev. Lett. 89, 195506 (2002). 

[20] G. Picard, A. Ajdari, F. Lequeux, and L. Bocquet, 
Phys. Rev. E 71, 010501 (2005). 























[21] M. Talamali, V. Petaja, D. Vandembroucq, and S.Roux, 
Phys. Rev. E 84, 016115 (2011). 

[22] J. Lin, A. Saade, E. Lerner, A. Rosso, and M. Wyart, 
Europhys. Lett. 105, 26003 (2014). 

[23] J. Lin, E. Lerner, A. Rosso, and M. Wyart, PNAS 111, 
14382 (2014). 

[24] J. Lin and M. Wyart, arXiv: 1506.03639 (2015). 

[25] G. Picard, A. Ajdari, F. Lequeux, and L. Bocquet, 
Eur. Phys. J. E 15, 371 (2004). 

[26] F. Puosi, J. Rottler, and J-L. Barrat, Phys. Rev. E 89, 
042302 (2014). 

[27] C. Maloney and A. Lemaitre, Phys. Rev. Lett. 93, 195501 
(2004). 

[28] A. Tanguy, B. Mantisi and M. Tsamados, Euro¬ 
phys. Lett. 90, 16004 (2010). 

[29] M. L. Manning and A. J. Liu, Phys. Rev. Lett. 107 
108302 (2011). 

[30] J. Rottler, S. S. Schoenholz, and A. J. Liu, Phys. Rev. E 
89, 042304 (2014). 

[31] S. Wijtmans, M. L. Manning, arXiv:1502.00685 (2015). 

[32] N. Xu, V. Vitelli, A. J. Liu, and S. R. Nagel, Euro¬ 
phys. Lett. 90 56001 (2010). 

[33] Finding its global minimum is desirable, but is currently 
a formidable task; we leave the design of algorithms that 
can find the global minimum of b(z) for future studies. 


[34] In a system of linear size L, the scaling law k,^, ~ y/8y is 
only expected to hold below 5y ~ L -4 due to hybridiza¬ 
tions with low-frequency plane waves. 

[35] C. S. O’Hern, L. E. Silbert, A. J. Liu, and S. R. Nagel, 
Phys. Rev. E 68, 011306 (2003). 

[36] A. J. Liu and S. R. Nagel, Annu. Rev. Condens. Matter 
Phys. 1, 347 (2010). 

[37] M. van Hecke, J. Phys. Condens. Matter 22, 033101 

( 2010 ). 

[38] We consider the absolute magnitudes |r^,| since normal 
modes T, and therefore also their associated asymmetries 
t$, are defined up to a sign. 

[39] E. Lerner, E. DeGiuli, G. During, and M. Wyart, Soft 
Matter 10, 5085 (2014). 

[40] Prestress is a generic feature of glassy solids; see, e.g., 
S. Alexander, Phys. Rep. 296, 65 (1998). 

[41] M. Wyart, L. E. Silbert, S. R. Nagel, and T. A. Witten, 
Phys. Rev. E 72, 051306 (2005). 

[42] http://www.netlib.org/lapack/ 

[43] http://mathworks.com/products/matlab/ 

[44] S. Karmakar, E. Lerner, and I. Procaccia, Phys. Rev. E 
82, 026105 (2010). 

[45] L. Gartner and E. Lerner, in preparation. 




